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1.  INTRODUCTION 

Sound  radiated  from  a  source  located  in  a  layer  of  fluid  of 
uniform  depth  overlying  a  second  fluid  with  greater  sound  speed  (a 
"fast"  bottom)  can  be  trapped  in  the  upper  layer.  This  trapping 
mechanism  can  be  described  by  either  ray  or  normal-mode  theory:  By 
ray  theory,  if  the  angle  of  incidence  of  the  ray  on  the  bottom 
exceeds  the  critical  angle  9C 

cos  0C  -  C1/C2 

where  9  is  the  grazing  angle  the  ray  makes  with  the  bottom  and 
ci  and  C2  are  the  speeds  of  sound  in  the  top  and  bottom  respec¬ 
tively,  the  acoustic  energy  is  totally  reflected  and  no  energy  is 
transmitted  into  the  bottom.  By  normal  mode  theory,  each  normal 
mode  has  associated  with  it  a  cutoff  frequency 

fn  *  (ci/2h) (n-l/2)cosec  9C 

where  n  is  the  mode  number  and  h  the  thickness  of  the  upper  layer. 
For  frequencies  above  cutoff,  the  mode  carries  energy  in  the  upper 
layer;  acoustic  energy  is  trapped  in  the  upper  layer  and  propagates 
to  great  distances.  If  both  boundaries  are  smooth,  the  propagation 
loss  is  only  that  associated  with  cylindrical  spreading  and 
absorption  losses. 

If  the  upper  layer  is  tapered  so  that  it  forms  a  wedge  subtend 
ing  an  angle  8,  then  as  sound  propagates  towards  the  shallower 
portion  of  the  layer  it  will  reach  a  depth  at  which  transmission  of 
acoustic  energy  into  the  lower  layer  is  possible.  From  the  ray 
point  of  view,  a  trapped  ray  will,  on  each  encounter  with  the 
sloping  bottom,  increase  its  angles  of  incidence,  until,  after  a 
sufficient  number  of  encounters,  the  grazing  angle  will  have 


increased  to  a  value  greater  than  critical  and  energy  enters  the 
bottom.  From  the  mod**  point  of  view,  as  the  depth  of  the  upper 
layer  decreases,  the  cutoff  frequency  of  a  given  mode  increases 
until  eventually  the  cutoff  frequency  exceeds  that  of  the  sound  and 
energy  enters  the  bottom.  For  a  given  frequency  f,  the  thickness  of 
the  wedge  at  which  the  lowest  mode  attains  cutoff  is 
h^  =  A/(4  sin  0C) 

where  \  -  cj/f  is  the  wavelength  in  the  wedge.  The  distance  from 
the  apex  measured  along  the  interface  at  which  the  lowest  mode 
attains  cutoff,  the  "dump  distance*  is 

xo  a  hj/tan3  *x/(4  sin  ©c  tanfj  )  (1) 

Obtaining  an  exact  solution  to  the  propagation  of  sound  from 
the  wedge  into  the  bottom  is  difficult  because  the  problem  is  not 
separable;  refraction  at  the  bottom  interface  results  in  components 
of  the  propagation  vector  at  the  bottom  that  depend  on  x.  This 
precludes  writing  the  boundary  condition  for  the  normal  component  of 
particle  velocity  at  the  interface  in  terms  of  one  coordinate.  As  a 
result,  the  acoustic  field  cannot  be  considered  as  a  product  of 
functions  of  separate  spatial  coordinates. 

Many  of  the  usual  approximate  techniques  also  become  inappli¬ 
cable.  For  example,  in  the  "adiabatic”  assumption  each  normal  mode 
in  the  wedge  is  assumed  to  retain  its  identity.  This  means  that  the 
problem  can  be  treated  using  normal  modes,  whose  propagation  vectors 
will  adjust  slowly  so  that  the  propagation  looks  as  if  it  were 
occurring  in  a  layer  of  constant  depth  whose  depth  is  the  local 
depth.  However,  when  cutoff  is  reached,  the  discrete  set  of  normal 
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modes  is  replaced  by  the  continuous  set  of  untrapped  modes,  and  the 
concept  of  modal  identity  loses  its  validity. 

Our  approach  is  to  use  the  method  of  images  to  predict  the 
amplitude  and  phase  of  the  pressure  at  the  bottom  interface.  The 
method  of  Green's  functions  is  then  used  to  obtain  the  beam  pattern 
of  the  sound  radiated  into  the  bottom.  The  reflections  of  sound 
from  the  surface  and  bottom  are  replaced  by  images  of  the  source  in 
the  same  way  as  for  a  sound  source  in  a  fluid  layer  overlying  an 
infinite  fluid  half-surface,  a  classical  problem  first  solved  by 
Pekeris.  In  the  case  of  a  wedge,  the  images  lie  on  the 
circumference  of  a  circle  whose  center  is  the  apex  of  the  wedge. 

The  lowest  images  (closest  to  the  source)  correspond  to  rays  of 
sound  which  make  grazing  reflections  from  the  surfaces  of  the  wedge. 
Higher  images  correspond  to  rays  with  greater  angles  of  elevation 
and  depression;  these  rays  suffer  more  reflections  from  the  surfaces 
of  the  wedge.  Finally,  images  are  encountered  for  which  the  rays 
exceed  the  critical  angle  at  the  bottom.  Still  higher  images 
correspond  to  more  reflections  from  the  bottom  (with  reflection 
coefficients  less  than  one  in  magnitude)  so  the  effective  strength 
of  these  higher  images  will  be  progressively  reduced.  (This  avoids 
certain  problems  involved  with  diffraction  from  the  apex  of  the 
wedge.)  Phase  coherent  summation  of  the  contributions  from  the 
source  and  images  yields  the  desired  pressure  and  phase  at  the 
bottom  of  the  wedge. 

Once  the  pressure  amplitude  and  phase  distribution  on  the 
interface  are  known,  the  interface  can  be  treated  as  an  apparent 
acoustical  source.  Application  of  straightforward  Green' s-f unction 
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techniques  then  results  in  an  integral  expression  for  the  sound 
field  in  the  bottom, 

^  '«  <J 


I A  -  A.  |  -  vW)  +  (§'$') 

Under  the  assumption  that  the  field  in  the  bottom  is  studied  for 
Ia.--a/I  X  ,  these  equations  can  be  well  approximated  by 

^  ,x  ITJX  ( TYf'.)  Jf' 

[(?. ry.fr 

$s/*A”  7sf/Xo 
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and  in  the  above  integrals 
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(2) 


,u>t  i[iot +Sl(f)  ] 
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P(J')e  -Prj')e 


?>o 


(3) 


is  the  image-calculated  pressure  at  the  interface  and  we  must  have 
P(£')  *  0  for  5'  <  0  to  satisfy  the  boundary  condition  to  the  right 
of  the  apex  of  the  wedge. 

The  first  stage  of  this  investigation,  the  subject  of  this 
report,  is  the  prediction  of  the  pressure  and  phase  distribution 
along  the  Interface. 
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THEORY 


a .  Finite  distance  from  source  to  receiver 
Figure  1  shows  a  wedge  with  a  pressure-release  surface  and  a 
fast  bottom.  The  origin  of  coordinates  is  at  the  apex  with  the  x 
axis  along  the  bottom,  positioned  at  the  left,  and  the  y  axis 
positioned  downward.  A  source  is  a  distance  x  from  the  axis  and  a 
depth  H  below  the  sea  surface.  The  line  connecting  the  source  to 
the  apex  makes  an  angle  <fd  with  the  sea  surface  and  an  angle 


(4) 


with  the  bottom.  The  distance  between  a  receiving  point  on  the 
bottom  at  distance  x  from  the  nfch  image  is 


*  \J  X  Z * *+ 1  ~  2/x.  X  c-o-s- 

(5) 

where  <^is  the  angle  between  the  line  joining  the  image  to  the  apex 
and  the  sea  bottom.  (See  Fig.  1.)  Each  is  given  by 

n 


where  INT[  ]  denotes  the  largest  integer  which  is  equal  to,  or 
smaller  than  the  argument.  The  number  of  times  an  image  interacts 
with  the  bottom  is  equal  to  the  number  of  times  the  line  from  the 
image  to  the  receiving  point  intersects  the  images  of  the  bottom. 
The  grazing  angle  on  the  bottom  for  the  mth  bounce  (counted  away 
from  the  receiving  point)  for  the  nfch  image  is  the  angle 
between  the  line  joining  the  nth  image  to  the  receiver  and  the 
m**1  image  of  the  bottom, 
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(7) 


$  . 9  -2m8 

r>  m  rtO  i 

where  Qno  is  the  grazing  angle  at  the  receiving  point. 

If  the  distance  from  the  source  to  the  point  on  the  bottom  at 
which  the  initial  bounce  takes  place  is  large  compared  to  a  wave¬ 
length  in  the  wedge,  the  reflection  coefficient  for  the  nfch 
image  at  the  mfc^  bounce,  can  be  approximated  by  the  plane-wave 
reflection  coefficient, 
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where  c^,  C2,  p\t  P2  a^e  the  sound  speeds  and  densities  of 
the  two  media,  and 
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For  a  source  of  unit  pressure  at  lm  from  its  center,  the 
complex  pressure  £(x,t)  =  P(x)e3ut  infinitesimally  above  the 
bottom  has  complex  amplitude 


P(V)  =  ±  ‘ 

•i  r 


r"  -i 

TT  R„„] 


+  l{iy  (-0  TT  5„„  ] 

hz  o  m  -<? 

where  the  first  term  is  the  contribution  from  the  direct  path,  the 
second  tt.-.m  is  the  contribution  from  the  images  above  the  wedge  and 
the  last  term  denotes  the  contribution  of  the  images  below  the 
wedge;  k  is  the  wave  number  in  the  wedge; 


M  , 

"  1  6 


r  TT  " 


II  -  1  N  I  - -  *  "  -  "  ‘  O  1 

n  1  {3  J  L  Y  (11) 

and  the  (-i)INTC  3  accounts  for  the  phase  shift  on  each  interaction 
with  the  pressure  release  surface.  The  second  and  third  terms 
differ  only  in  that  there  is  a  0^  bounce  associated  with  the 
images  above  the  wedge. 

The  complex  pressure  amplitude  along  the  bottom  can  be  written 


as  a  simple  sum 


pw  =  L£e 

— ■*  V) 


)TT  T? 

—  no 


Computer  program  (WEDGE  1)  calculates  the  pressure-amplitude  and 
phase-distribution  along  the  bottom  resulting  from  a  source  at  a 
finite  distance  from  the  receiver. 
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A  refinement  which  includes  the  effect  of  a  directional  source 


is  presented  in  M.  Kawamura  and  I.  Ioannou,  "Pressure  on  the  Inter¬ 
face  Between  a  Converging  Fluid  Wedge  and  a  Fast  Fluid  Bottom",  M.S. 
Thesis,  NPS ,  (Dec.  1978)  and  named  WEDGE  0.  In  both  of  these  pro¬ 
grams  the  complex  pressure  amplitude  P(x)  was  "normalized"  by  multi¬ 
plying  the  right  hand  side  of  Eq(9)  by  rQ.  This  has  the  effect  of 
assuring  that  the  source  generates  a  pressure  amplitude  at  one  meter 
which  is  proportional  to  its  distance  from  the  apex.  The  advantage 
is  that  the  effects  on  rQp(x)  of  changing  rQ  are  easily  studied 
without  having  to  correct  for  the  l/rQ  diminuation  in  P(x). 
b.  Infinite  distance  from  source  to  receiver. 

In  the  case  of  a  distant  source,  inspection  of  Fig.  2  yields  the 
approximations 


9*. i  4>„ 

Now,  since 


(13) 


and  the  expressions  for  the  reflection  coefficients  are 

TT  8... =  f „  *  B  (O 

v*  -o 

TT  1?,  *1?,  -RC4!) 

fn?  ,JT  R„  1R  *  TT£(<P  )  R(4>J 


(14) 


(15) 


X 


Fig.  2  Geometry  of  the  distant 


The  right  hand  sides  of  Eq.  15  depend  only  on  n,  so  we  can  define 


r  m 


!  Z. 

)  B. 

r  (—1')  R, 


n  =  0 

*  =  f 

*  *  Z,  V  #  * 


Also,  the  geometry  of  Fig.  2  gives 

rn  =  r  -  x  cos  <fn 

Substitution  of  Eqs .  (15),  (16)  and  (17)  into  Eq .  (12)  yields 


(16) 


(17) 


(18) 
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Computer  program  (WEDGE  2)  evaluates  Eq.  (18)  and  can  be  used  to 
calculate  the  pressure  amplitude  and  phase  distribution  along  the 
bottom  resulting  from  a  source  infinitely  distant  from  the 
receiver. 

As  in  WEDGE  1,  what  was  actually  computed  was  the  "normalized" 
complex  pressure  amplitude  rP(x) . 

3.  TESTING. 

a .  Verification  of  WEDGE  1 

To  verify  that  WEDGE  1  was  performing  the  calculations  correct¬ 
ly,  several  runs  were  made  for  bottoms  with  high  impedances  to 
approximate  a  rigid  boundary.  For  these  runs  the  wedge  -angle  was 
chosen  to  produce  a  limited  number  of  images  so  that  calculation  of 
the  pressure  was  simple  enough  to  be  performed  on  a  hand  calculator. 

Runs  were  made  for  90°  and  45®  wedge  angles.  In  all  cases  the 
agreement  between  hand  calculation  and  WEDGE  1  output  was  considered 
to  be  very  good.  An  example  is  given  in  Figs.  3  and  4  for  a  45s 
wedge  Cj/C2  *  0.2  and  p./p,.  *  0.05.  In  these  and  the  following 
figures  the  pressure  amplitude  is  normalized  so  that  its  greatest 
value  is  unity.  (This  normalization  was  not  done  in  WEDGE  0.) 

b.  Verification  of  WEDGE  2 

Having  affirmed  that  WEDGE  1  was  producing  valid  output  for  the 
pressure  amplitude  and  phase  distribution  at  the  interface  between 
the  wedge  and  the  bottom,  comparisons  between  WEDGE  1  and  WEDGE  2 
were  made.  Because  WEDGE  2  is  valid  only  for  distant  sources, 
comparisons  were  performed  for  source  distances  greater  than  100 
times  the  dump  distance  X0.  One  example  of  this  comparison  is 
given  in  Fig.  5.  WEDGE  2  results  are  plotted  only  where,  by  visual 
inspection,  they  do  not  exactly  overlap  the  results  of  WEDGE  1. 
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Normalized  Distance  (x/X^) 
Fig.  4  Comparison  of  phase  angles 
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Fig.  5  Comparison  of  pressure  amplitudes 


For  the  WEDGE  1  calculations  the  source  was  located  647  dump 
distances  from  the  apex  and  at  mid-depth.  For  WEDGE  2,  the  source, 
defined  to  be  at  infinity,  was  positioned  at  mid-depth  so  that  the 
same  normal  modes  were  excited  in  each  case. 

In  all  cases  investigated,  the  amplitudes  of  the  pressure  on 
the  bottom  calculated  by  the  two  programs  were  in  excellent 
agreement. 

Phase  results  are  not  displayed  because  the  phases  fluctuated 
so  rapidly  with  respect  to  the  increment  in  distance  that  no 
discernable  pattern  could  be  seen.  However,  detailed  comparisons  of 
the  numerical  outputs  of  WT'OGE  1  and  WEDGE  2  confirmed  that  the  two 
programs  were  also  in  agreement  in  their  calculations  of  the  phases 
of  the  pressures  at  >_he  interface. 

4 .  DISCUSSION  OF  RF/l  H  fjj 

Figures  5  and  6  show  the  pressure  amplitude  (normalized  to  give 
a  peak  value  of  unity)  as  a  function  of  distance  from  the  apex 
(normalized  by  dividing  by  the  dump  distance) .  Both  figures  are  for 
the  same  wedge  geometry  and  wedge  and  bottom  properties?  they  differ 
only  in  the  location  of  the  source.  In  Fig.  5  the  WEDGE  1,  source 
is  at  647  dump  distances  and  at  mid-depth.  In  Fig.  6  WEDGE  1,  the 
source  is  at  64.7  dump  distances  and  at  a  depth  of  1000  m  (one-fifth 
of  the  depth  at  that  distance) .  In  both  figures  the  WEDGE  2 
calculation  was  carried  out  with  the  source  position  to  give  the 
same  source  angle  as  the  corresponding  WEDGE  1  calculations. 


16 


In  all  cases  the  pressure  amplitude  on  the  bottom  increases 
smoothly  from  zero  at  the  apex  to  a  local  maximum  at  the  dump 
distance.  The  pressure  amplitude  then  decreases  at  distances 
further  from  the  apex  until  modal  interference  becomes  important  at 
about  two  dump  distances  when  the  position  of  the  source  is  such 
that  the  second  mode  is  strongly  excited  (Fig.  6)  and  at  about  three 
dump  distances  when  the  second  mode  is  weakly  excited  (Fig.  5). 

Comparison  of  the  results  for  WEDGE  1  and  WEDGE  2  in  Fig.  6 
shows  that  there  is  little  difference  between  the  pressure 
distribution  produced  by  a  source  at  65  dump  distances  and  one 
infinitely  distant;  the  positions  and  magnitudes  of  the  the 
interferences  peaks  are  slightly  different,  but  the  results  out  to 
about  two  dump  distances  are  indistinguishable  in  the  two  cases. 

5.  CONCLUSIONS 

Two  programs  have  been  generated  to  calculate  the  pressure 
amplitude  and  the  phase  distribution  on  the  interface  between  a 
fluid  wedge  and  an  underlying  semi-infinite  fluid  medium.  The  first 
program,  WEDGE  1,  is  the  more  general  in  that  it  allows  the  source 
to  be  positioned  anywhere  within  the  wedge.  The  second  WEDGE  2  is 
more  restricted  in  that  the  source  is  assumed  to  be  infinitely 
distant  from  the  apex  of  the  wedge.  In  practice  this  means  that  for 
reasonably  good  predictions  of  the  pressure  and  phase  at  the  inter¬ 
face  out  to  a  couple  of  times  the  dump  distance  XQ,  the  source 
should  lie  more  than  about  100  XQ  away  from  the  apex.  In  running 
comparable  cases,  the  difference  in  computation  time  between  the  two 
programs  amounts  to  a  factor  of  five,  with  WEDGE  2  being  the  faster 
of  the  two. 
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APPENDIX  I 


Equations  for  the  computer  program  (wedge  1) 


n  =  4  tan  ( 1) 

K  =  2  n  f  /  G 

Y  =  p  C 

9o=  fi  -  sin  (H/X) 

Pn  =  2  INT  (  (n+l)/2  ]/J  +  (-1)" 
$kT  tan'lsin^n/  (cos  <Pn  -  x/x)) 


constant 
wave  number 
acoustic  impedance 


Dna  'I  X2+  x2-  2  x  X  cos^n/x 
M  =  MT  (^n/2/5) 

(tom*  ftio”  2  m  /5 

=  V1  "  (P2/C1)  costtnm/  slndnm 

_ 

^2  =  -  i  f(C2/C1)2cos ^  1  ~  /  slndnr 

or  & 

5nf  (%/r,  *  <jy  /  (r2/r,  ♦  ^ 

P  *  ^1  1/Bn-  .-il'xD"(-l)'N,'^L’  ( I 

P-I£l 

/P  =  tan’l  Im(P)/Re(P)  |  pressure  phase 


angle  between  source  and 
bottom 

angle  between  source  and 
n'th  image 

incident  angle  to  the  bottom 
(last  reflection) 
normalized  distance  between 
image  and  receiver 
bottom  reflections 


incident  angle  to  the  bottom 
(intermediate  reflections) 
parameter  for  a  reflection 
coefficent  (  0rrr^  %) 
parameter  for  a  reflection 
coefficent  (e^  6>c) 


reflection  coefficent 


mo  pressure  at  point  (x,0) 
pressure  amplitude 


*  INT( )  max.  integer  of  an  argument 

2  INT  (n+l)/2  =  n  +  l/2[l  -  (-1)"] 

**  INT  (n+l)/2  =  1/4  [2n  +  1  -  (-1  f  } 
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APPENDIX  II 
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C  WEDGE  1 
C  ***  MAIN  PROGRAM  *** 

C 

C  ********************************************************************* 

C  *  THIS  PROGRAM  EMPLOYS  THE  METHOD  OF  IMAGES  TO  OBTAIN  THE  PRESSURE  * 
C  *AND  PHASE  DISTRIBUTION  ALONG  THE  BOTTOM  OF  A  WEDGE-SHAPED  FLUID  * 


C  *LAYER. OVERLYING  A  FAST  BOTTOM.  * 
C  *  * 
C  *  ASSUMPTIONS  * 
C  *  (1)  PLANE  WAVE  REFLECTION  COEFICIENTS.  * 
C  *  (2)  SOURCE  IS  POINT  SOURCE.  * 


C  ********************************************************************* 

c 

C  *  NOTATION 


c 

I/O 

OUT-OUTPUT 

c 

IN  “INPUT 

c 

CHARACTER 

c 

V  “VARIABLE 

c 

P  “PARAMETER 

c 

C  “CONSTANT 

c 

R  -REAL  NUMBER 

c 

I  “INTEGER 

c 

c 

c 

Z  “COMPLEX  NUMBER 

c 

c 

SYMBOL 

MEANING 

I/O 

CHARACTER 

c 

PA 

PRESSURE  AMPLITUDE  ON  THE  BOTTOM 

OUT 

V.R 

c 

PH 

PRESSURE  PHASE  ON  THE  BOTTOM 

OUT 

VrR 

c 

Cl 

SOUND  SPEED  IN  MEDIUM  1  (WEDGE) 

IN 

C  t  R 

c 

C2 

SOUND  SPEED  IN  MEDIUM  2  (BOTTOM) 

IN 

CrR 

c 

ROl 

DENSITY  OF  MEDIUM  1 

IN 

C.R 

c 

R02 

DENSITY  OF  MEDIUM  2 

IN 

C.R 

c 

BETA 

SLOPE  ANGLE  OF  BOTTOM 

IN 

CrR 

c 

X 

SOURCE  DISTANCE  FROM  THE  APEX  IN  METERS 

IN 

CrR 

c 

XX 

DISTANCE  FROM  THE  APEX  IN  METER 

OUT 

VrR 

c 

AT  UHICH  PA  AND  PH  ARE  CALCULATED 

c 

F 

DRIVING  FREQUENCY 

IN 

CrR 

c 

H 

SOURCE  DEPTH  IN  METERS 

IN 

C.R 

c 

ANGLC 

CRITICAL  ANGLE 

CrR 

c 

HO 

LOWEST  POSSIBLE  MODE  CUT  OFF  DEPTH  IN  METERS 

C.R 

c 

XO 

LOWEST  POSSIBLE  MODE  CUT  OFF  DISTANCE  FROM  THE 

C.R 

c 

APEX  IN  METERS 

c 

XI 

DISTANCE  FROM  APEX  AT  WHICH  CALCULATION  STARTS 

IN 

CrR 

c 

X2 

DISTANCE  FROM  APEX  AT  WHICH  CALCULATION  STOPS 

IN 

CrR 

c 

DX 

RANGE  OF  CALCULATION  IN  METERS 

C.R 

c 

N 

NUMBER  OF  POINTS  FOR  WHICH  PA  AND  PH  ARE 

IN 

C.I 

c 

CALCULATED 

c 

PA  I 

3.1415 . . 

CrR 

c 

WN 

WAVE  NUMBER 

CrR 

c 

C21 

SOUND  SPEED  RATIO  “C2/C1 

CrR 

c 

RC21 

ACOUSTIC  IMPEDANCE  RATIO 

C.R 

c 

ANGLO 

ANGLE  FORMED  BY  SOURCE r  APEX  AND  BOTTOM 

CrR 

c 

DISTX 

TEMPORARY  FOR  XX(I)/X 

PrR 

c 

PAMP 

TEMPORARY  FOR  PA 

PrR 

c 

PHAS 

TEMPORARY  FOR  PH 

PrR 

c 

C 

C 


C  ********************************************************************** 
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oonn  oono  ooonoo  nno  o  ci  o  nn non 


( (PROGRAM) ) 


DEFINE  CHARACTERS  OF  VARIABLES. PARAMETERS  AND  CONSTANTS 

IMPLICIT  REAL*8(A-H.0-Z) . INTEGER ( I-N) 

DIMENSION  PA< 101 ) »PH< 101 ) . XX < 101) 

COMMON  C21 » RC2 1 . BETA . UN  >  PA I . ANGL  0 • X  >  H • F 

READ  IN  DATA 

READ(S.IOO)  C1.C2.R01.R02. BETA 
READ <3. 500)  X1.X2.X.H.F 
READ (5. 600)  N 

CALCULATE  AND  URITE  CONSTANTS 

WRITE  (6.700)  Cl  .C2.R01  >R02.  XI  .X2. X.H. BETA.F 
ANGLC-DARCOS ( C1/C2 ) 

H0-C1/(4.0D0*F*DSIN( ANGLO ) 

XO-HO/DSIN( BETA ) 

WRITE (6 » 400 )  HO.XO 

CALCULATE  RANGE  DX 

DX-X2-X1 

CALCULATE  CONSTANTS  FOR  SUBROUTINE 

PAI-4 . 0D0*D AT AN ( 1 . 0D0 ) 

UN»2.0D0*PAI*F/C1 

C21-C2/C1 

RC21-C2*R02/(C1*R01) 

ANGLO-BETA-DARSIN ( H/X ) 

URITE  (6.200) 

CALCULATE  PRESSURES  WITH  RESPECT  TO  DISTANCE  XX 
WHERE  XX  IS  CHANGED  BY  INCRIMENT 

L-N+l 

DO  101*1 .L 

XX ( I ) -XI +DFLOAT ( I - 1 >  *DX/DFLOAT ( N ) 

DISTX-XX ( I ) /X 

CALL  PRES  (DISTX.PAMP.PHAS) 

PA  < I ) -PAMP 
PH ( I ) -PHAS 

URITE  (6.300)  XX ( I ) »PA( I ) »PH( I > 

10  CONTINUE 

URITE  (6.800) 

PLOT  THE  PRESSURE  AMPLITUDE  VS.  DISTANCE  XX 
AND  PRESSURE  PHASE  VS.  DISTANCE  XX 

URITE  (6.700)  Cl. C2.R01.R02. XI. X2. X.H. BETA.F 
CALL  DPLTP  (XX.PA.N.O) 

URITE  (6.800) 

URITE  (6.700)  Cl .C2.R01 .R02.X1 .X2. X.H. BETA.F 
CALL  DPLTP  (XX.PH.N.O) 
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C  I/O  FORMATS 
C 

100  FORMAT  (5F13.0) 

200  FORMAT  <//22X. 'DISTANCE' » 17X. 'PRESSURE  AMP. '» 12X» 'PHASE  ANGLE'//) 
300  FORMAT  (  20X.D15. 7. 10X,015. 7. 10X. D1S. 7) 

400  FORMAT  {//20X» 'LOWEST  POSSIBLE  MODE  CUT  OFF  DEPTH' »9X>D15.7 
C  /20X. 'DISTANCE  FROM  APEX  ' »20X*D15. 7//) 

500  “FORMAT  TSFl  5  rO) 

600  FORMAT  (13) 

700  FORMAT  <//25X. ' INPUT  DATA' . 13X. ' UEDGE1 '// 

C  25X.'C1  *' .D15.7.3X. 'C2  -'»D15.7/ 

C  25X.'R01  *' .D15.7.5X. 'R02  -'.DlS.7/ 

C  25X.'X1  *' .D13.7.5X. 'X2  -'*D15.7/ 

C  25X.'X  -'  .D15.7.5X.  'H  »'»D15.7/ 

C  25X»'B£TA-'.D13.7»5X.'F  -'»D15.7/> 

800  FORMAT  <1H1> 

STOP 

END 
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nnnnnnnnnnnnnnoftoonnnnftnnnnnnonnnnnnnnnn 


***  SUB  PROGRAM  *** 

********************************************************************* 

*  SUBROUTINE  'PRES'  CALCULATES  THE  PRESSURE  AMPLITUDE  AND  PHASE  AT  A* 

•POINT  ON  THE  BOTTOM  BY  USING  METHOD  OF  IMAGES.  * 

********************************************************************* 

UHERE  C21 .RC21 .  BETA .UN .PA I 1X1H  AND  F  ARE  INPUT  CONSTANTS 
FROM  MAIN  ROUTINE. 

DISTX  IS  INPUT  VARIABLE  FROM  MAIN  ROUTINE.  AND. 

PAMP  AND  PHAS  ARE  OUTPUTS  TO  MAIN  ROUTINE. 

*  NOTATION 


SYMBOL  MEANING  I/O  CHARACTER 


PRES  PRESSURE  CAUSED  BY  N  TH  PAIR  OF  IMAGES  P.Z 

P  SUM  OF  PRES. (IF  ALL  POSSIBLE  PAIR  OF  IMAGES  ARE  P.Z 

ADDED  P  IS  PAMP 

ANGLN  ANGLE  FORMED  BY  IMAGE.  APEX  AND  BOTTOM  P.R 

DD  PARAMETER  FOR  CONVENIENCE  P.R 

R  TOTAL  REFLECTION  LOSS  OF  <N'>  TH  IMAGE  P.Z 

N  N  TH  IMAGE  OR  N  TH  PAIR  OF  IMAGES  PrI 

M  M  TH  BOTTOM  BOUNCE  P.I 

MM  MAXIMUM  NUMBERS  OF  BOTTOM  BOUNCES  P.I 

THETAO  ANGLE  FORMED  BY  THE  N  TH  IMAGE. THE  RECEIVING  POINT 

AND  THE  BOTTM  P.R 

DR  DISTANCE  FROM  N  TH  IMAGE  TO  THE  RECEIVING  POINT 

NORMALIZED  BY  DIVIDING  BY  XO  P.R 

THETAM  GRAZING  ANGLE  TO  THE  BOTTOM  P.R 

(INTERMEDIATE  REFLECTION) 

CHECK  IDENTIFIES  IF  GRAZING  ANGLE  EXCEEDS  CRITICAL  ANGLE  P.R 

PSA I  PARAMETER  IN  REFLECTION  COEFFICIENT  P.Z 

REFL  REFLECTION  COEFFICIENT  OF  INTERMEDIATE  BOUNCE  P.Z 

REFLNO  REFLECTION  COEFFICIENT  OF  THE  LAST  BOUNCE  P.Z 

2  PARAMETER  FOR  CONVENIENCE  P.Z 


********************************************************************** 


nnnnn  onoo  non  ono  non  oon  noon  non  ono  non  on 


( (PROGRAM) ) 

SUBROUTINE  PRES  (DISTX.PAHP.PHAS) 

DEFINE  CHARACTERS  OF  VARIABLE t PARAMETERS  AND  CONSTANTS 

IMPLICIT  REAL#8(A-H»0-Z> » INTEGER < I-N> 

COMPLEX* 1 6  PSA  I »ft£FL> REPLNO , R , Z t PRES . P 
COMMON  C21 »  RC21 »  BETA . UN  t P A I , ANGLO . X » H  »  F 

INITIALIZE  COMPLEX  PRESSURE  P-O.O+JO.O 

P«DCMPLX<O.ODO»O.ODO> 


RESET  COUNTER  N-0 
N-0 

10  CONTINUE 

CALCULATE  ANGLN  AND  DETERMINE  IP  ANGLE  IS  IN  THE  RANGE  OR  NOT 
IP  ANGLE. GT.PAI  GO  TO  THE  NEXT  STEP 

ANGLN-2.0D0*INT  < (N+t . >/2. >*BETA+(-l )**N*ANGL0 
IF  < ANGLN. GE. PAI)  GO  TO  SO 

CALCULATE  PARAMETER  THETAO  AND  DR 

DD-DCOS< ANGLN )-DISTX 
THETA0-DATAN2<DSIN(ANGLN> »DD> 

DR-DSQRT ( 1 . ODO+D ISTX**2-2 . ODO*DISTX*DCOS ( ANGLN >  > 

INITIALIZE  THE  TOTAL  REFLECTION  COEFFCIENT  R-1.0+J0.0 
R-DCMPLX< 1 .ODO.O.ODO) 

DEFINE  MM  BY  EACH  PATH  UHICH  IS  DETERMINED  BY  THE  N  TH  IMAGE 

MM-IDINT  <  ANGLN/ <  2  <  0D04BETA ) ) 

L-MM+1 
00  40  M-l  »L 

CALCULATE  PARAMETERS  THETAM  AND  CHECK 

THETAM-THETAO-2 . 000*DFL0AT<M-1 >  *BETA 
CHECK- 1 . 0D0-C21 **2*DC0S <  THETAM  >  **2 

IDENTIFY  THE  INCIDENT  ANGLE  UHICH  IS  LESS  THAN  THE  CRITICAL 
ANGLE  OR  NOT 

IF<CHECK.GT. 0.000)  GO  TO  20 

CALCULATE  PARAMETER  PSAI  FOR  THE  REFLECTION  COEFFCIENT  OF  EACH 
BOUNCE.  THERE  ARE  THO  WAYS  UHICH  DEPEND  ON  THE  IDENTIFICATION 
OF  A  CHECK 

PSA I-DCMPLX  <  0 . 000 t -DSORT  < -CHECK ) /DS I N  <  THETAM ) ) 

GO  TO  30 
20  CONTINUE 

PSAI-OCMPLX<OSORT(CHECK)/DSIN<  THETAM) » 0.000) 

30  CONTINUE 
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-  -A.  ^ . 


C  CALCULATE  REEL 
C 

REFL • <  RC2 1 -PS A I ) / <  RC2 1 +PS A I ) 

C 

C  IF  THE  BOUNCE  IS  THE  LAST  FROM  THE  SOURCE » THEN  REFLNO-REFL 
C 

IF(M.EQ. 1 >  REFLNO-REFL 
C 

C  CALCULATE  THE  TOTAL  REFLECTION  COEFFCIENT  R  BY  EACH  PATH 
C 

R-RBREFL 
40  CONTINUE 
C 

C  CALCULATE  PARAMETERS  Z  AND  PRES 
C 

Z-DCMPLX <  0 . ODO . -WN*X*DR  > 

PRES-CDEXP  < ZTfr=l  Tf*THTC< hffl  7 >72. )  *' 

C  <l.ODO+1.0DO/REFLNO)*R/DR 

C 

C  CALCULATE  THE  PRESSURE  P  WHICH  IS  A  SUH  OF  PRES 
C 

P-P+PRES 

C 

C  SET  TO  THE  NEXT  PAIR  OF  IMAGE 
C 

N-N+l 
00  TO  10 
SO  CONTINUE 
C 

C  CALCULATE  PA  AND  PH. AND  RETURN  BACK  TO  THE  MAIN  ROUTINE 
C 

PAMP-CDABS ( P ) 

PREAL-P 
P— P— PREAL 

PP-PBDCMPLX ( 0 . ODO . 1 . ODO  > 

PIMAG-DSIGN<CDABS(P) .-PP) 

PHAS-DAT AN2  <  P I MAG . PREAL ) 

RETURN 

END 


(WEDGE  2) 

• PRKS  * 
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ooonnnonnoo nnn  nonooooooooooonooooooooooonnnooooooooooonooooo 


WEDGE  2 

*»*  MAIN  PROGRAM  *** 

******************************************************************** 

*  THIS  PROGRAM  EMPLOYS  THE  METHOD  OF  IMAGES  TO  OBTAIN  THE  PRESSURE 
*AND  PHASE  DISTRIBUTION  ALONG  THE  BOTTOM  OF  A  UEDGE-SHAPED  FLUID 
♦LATER » OVERLYING  A  FAST  BOTTOM. 

* 

*  ASSUMPTIONS  > 

*  <1>  PLANE  WAVE  REFLECTION  CQEFICIENTS. 

*  <2>  SOURCE  IS  A  POINT  SOURCE.  > 

*  (3)  SOURCE  IS  A  DISTANT  SOURCE. 

********************************************************************: 

*  NOTATION 

I/O  OUT-OUTPUT 
IN  -INPUT 
CHARACTER 

V  -VARIABLE 
P  -PARAMETER 
C  -CONSTANT 

R  -REAL  NUMBER 
I  -INTEGER 
Z  -COMPLEX  NUMBER 


SYMBOL 

MEANING 

M 

X 

Ol 

CHARACTER 

PA 

PRESSURE  AMPLITUDE  ON  THE  BOTTOM 

OUT 

VrR 

PH 

PRESSURE  PHASE  ON  THE  BOTTOM 

OUT 

VfR 

Cl 

SOUND  SPEED  IN  MEDIUM  1  (WEDGE) 

IN 

C.R 

C2 

SOUND  SPEED  IN  MEDIUM  2  (BOTTOM) 

IN 

C»R 

ROl 

DENSITY  OF  MEDIUM  1 

IN 

CfR 

R  02 

DENSITY  OF  MEDIUM  2 

IN 

CfR 

BETA 

SLOPE  ANGLE  OF  BOTTOM 

IN 

CfR 

XX 

DISTANCE  FROM  THE  APEX  IN  METER 

AT  WHICH  PA  AND  PH  ARE  CALCULATED 

OUT 

V.R 

F 

DRIVING  FREGUENCY 

IN 

C»R 

ANGLC 

CRITICAL  ANGLE 

CfR 

HO 

LOWEST  POSSIBLE  MODE  CUT  OFF  DEPTH  IN  METERS 

CfR 

XO 

LOWEST  POSSIBLE  MODE  CUT  OFF  DISTANCE  FROM  THE 
APEX  IN  METERS 

CfR 

XI 

DISTANCE  FROH  APEX  AT  WHICH  CALCULATION  STARTS 

IN 

"CTfir 

X2 

DISTANCE  FROM  APEX  AT  WHICH  CALCULATION  STOPS 

IN 

CfR 

DX 

RANGE  OF  CALCULATION  IN  METERS 

CfR 

N 

NUMBER  OF  POINTS  FOR  WHICH  PA  AND  PH  ARE 
CALCULATED 

IN 

CfI 

PAI 

3.1413 . 

CfR 

WN 

WAVE  NUMBER 

CfR 

C21 

SOUND  SPEED  RATIO  -C2/C1 

CfR 

RC21 

ACOUSTIC  IMPEDANCE  RATIO 

CfR 

ANGLO 

ANGLE  FORMED  BY  SOURCE f  APEX  AND  BOTTOM 

IN 

CfR 

DISTX 

TEMPORARY  FOR  XX(I)/X 

PfR 

PA  MP 

TEMPORARY  FOR  PA 

PfR 

PHAS 

TEMPORARY  FOR  PH 

PfR 

********************************************************************** 
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noon  oono  ono  ooo  ooo  o  n  n  oooo 


((PROGRAM)) 

DEFINE  CHARACTERS  OF  VARIABLES .PARAMETERS  AND  CONSTANTS 

IMPLICIT  REAL*8(A-H.0-Z) , INTEQER( I-N) 

DIMENSION  PA(101).PH(101),XX(101> 

COMMON  C2 1 * RC21 * BETA  * UN  * PAI * ANGL 0 . F 

READ  IN  DATA 

READ(5'100)  C1*C2*R01.R02*BETA 
READ(5*500>  XI .X2.F. ANGLO 
READ(5.600)  N 

CALCULATE  AND  WRITE  CONSTANTS 

WRITE  (6*700)  Cl .C2.R01 >R02*  XI *X2*  BETA.F * ANGLO 
ANGLC-DARCOS ( C1/C2 ) 

H0»C1/(4.0D0*F*DSIN(ANGLC) ) 

X0*H0/ DSIN(BETA) 

WRITE (6*400)  H0*X0 

CALCULATE  RANGE  DX 

DX-X2-X1 

CALCULATE  CONSTANTS  FOR  SUBROUTINE 

PAI-4 . 0D04DAT AN ( 1 . ODO ) 

WN»2.0D0*PAI*F/C1 

C21-C2/C1 

RC21-C2*R02/(C1*R01 ) 

WRITE  (6*200) 

CALCULATE  PRESSURES  WITH  RESPECT  TO  DISTANCE  XX 
WHERE  XX  IS  CHANGED  BY  INCRIMENT 

L-N+l 

DO  10  I-1,L 

DISTX-X1+DFL0AT( 1-1 )*DX/DFL0AT(N) 

XX(I)«DISTX 

CALL  PRES  ( DISTX  >  PAMP  *PHAS  > 

PA(I)*PAMP 
PH( I ) *PHAS 

WRITE  (6.300)  XX(I)*PA(I)*PH(I> 

10  CONTINUE 

PLOT  THE  PRESSURE  AMPLITUDE  VS.  DISTANCE  XX 
ANir  PRESSURE  PHASE^S.'  DISTANCE'XX 

WRITE  (6*800) 

WRITE  (6*700)  Cl >C2*R01 *R02> XI *X2* BETA* F> ANGLO 
CALL  DPLTP  (XX.PA.N.O) 

WRITE  (6*800) 

WRITE  (6*700)  Cl *C2 *R01 *R02 *X1 * X2* BETA >F* ANGLO 
CALL  DPLTP  (XX. PH. N*0) 
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on 


I/O  FORMATS 


100  FORMAT 
200  FORMAT 
300  FORMAT 
400  FORMAT 
C 

500  FORMAT 
400  FORMAT 
700  FORMAT 
C 
C 

c 

c 

c 

800  FORMAT 
STOP 


<//22X?*  DISTANCE' » 17X» 'PRESSURE  AMP. ' >12X» 'PHASE  ANGLE'//) 
<  20X.D15.7.  10X.D15.7*  10X.D15. 7)  „„_u<  _  - 

<//20X. 'LOWEST  POSSIBLE  MODE  CUT  OFF  DEPTH  .9X.D15.7 
/20X. 'DISTANCE  FROM  APEX  ' .20X.015.7//) 

(4F15.0) 

<//15X» ' INPUT  DATA' »25X. 'UEDQE2'// 

25X. 'Cl  *'.01S.7.5X»'C2  *'.015.7/ 

25X.'R0I  ■'.D15.7.5X. 'R02  *'»D15.7/ 

23X»'X1  »' .D15.7  »5X. ' X2  -'.D15.7/ 

25X. 'BETA*' .D15 .7.5X. 'F  *' .015.7/ 

42X.' SOURCE  ANGLE-' .015.7/ > 

( 1H1 ) 


C  ***  SUB  PROGRAM  **t 
C 

C  ********************************************************************* 
C  *  SUBROUTINE  'PRES'  CALCULATES  THE  PRESSURE  AMPLITUDE  AMD  PHASE  AT  A* 
C  ftPOINT  ON  THE  BOTTOM  BY  USING  METHOD  OF  IMAGES.  * 

C  a******************************************************************** 

c 

C  UHERE  C21.RC21. BETA. WN.PAI. ANGLO  AND  F  ARE  INPUT  CONSTANTS 
C  FROM  MAIN  ROUTINE. 

C  DISTX  IS  INPUT  VARIABLE  FROM  MAIN  ROUTINE.  AND. 

C  PAMP  AND  PHAS  ARE  OUTPUTS  TO  MAIN  ROUTINE. 

C 

C  *  NOTATION 
C 


SYMBOL  MEANING  I/O  CHARACTER 


PRES  PRESSURE  CAUSED  BY  N  TH  PAIR  OF  IMAGES  P.Z 

P  SUM  OF  PRES. (IF  ALL  POSSIBLE  PAIR  OF  IMAGES  ARE  P.Z 

ADDED  P  IS  PAMP 

ANGLN  ANGLE  FORMED  BY  IMAGE.  APEX  AND  BOTTOM  P.R 

N  N  TH  IMAGE  OR  N  TH  PAIR  OF  IMAGES  P.I 

I FLAG  IFLAG-0  WHEN  N  IS  EVEN 

IFLAG-1  WHEN  N  IS  ODD  P.I 

REFL  REFLECTION  COEFFICIENT  OF  INTERMEDIATE  BOUNCE  P.Z 

FUNCTION  PROGRAM  REFL  (ANGLN) 

CHECK  PARAMETER  FOR  FUNCTION  REFL  P.R 

(IDENTIFY  IF  GRAZING  ANGLE  EXCEEDS  CRITICAL  ANGLE) 

PSA I  PARAMETER  FOR  FUNCTION  REFL  P.Z 

REFL NO  REFLECTION  COEFFICIENT  OF  THE  LAST  BOUNCE  P.Z 

R  TOTAL  REFLECTION  LOSS  OF  (N'>  TH  IMAGE  P.Z 

RA  TEMPORARY  FOR  R  WHEN  N  IS  EVEN  P.Z 

RB  TEMPORARY  FOR  R  WHEN  N  IS  ODD  P.Z 

Z  PARAMETER  FOR  CONVENIENCE  P.Z 


C 

U  *************¥******************************»«***********«************ 
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non  oono  non  nnn  nnn  non 


< (PROGRAM) ) 

SUBROUTINE  PRES  (DISTX.PAMP.PHAS) 

DEFINE  CHARACTERS  OF  VARIABLE » PARAMETERS  AND  CONSTANTS 

IMPLICIT  R£AL*8< A-H»0-Z) » INTEOER( I-N) 

COMPLEX* 1 A  R » RA , RB » REFL  » REFLNO . P . PRES . Z 
COMMON  C21 .RC21 »8ETA.WN*PAI » ANGLOrF 

INITIALIZE  COMPLEX  PRESSURE  P-O.O+JO.O 

P-DCMPLX (O.ODO.O.ODO) 

RESET  COUNTER  N-0  ►  RESET  IFLAG-0 

N«0 

I FLAG*0 
10  CONTINUE 

CALCULATE  AN6LN  AND  DETERMINE  IF  ANGLE  IS  IN  THE  RANGE  OR  NOT 
IF  ANGLE. GT.PAI  GO  TO  THE  NEXT  STEP 

ANGLN*2 . 000* I  NT  < <N+1. >/2. ) *B£TA+<-l ) **N*ANGLO 
IF  (ANGLN.GE.PAI)  GO  TO  80 

CALCULATE  THE  TOTAL  REFLECTION  COEFFICIENT  R  BY  EACH  PATH 

IF  (IFLAG.EQ.l)  GO  TO  40 
REFLNO»REFL ( ANGLN  > 

IF  <N.GE.2>  GO  TO  20 

RA-REFLNO 

R«RA 

GO  TO  30 
20  CONTINUE 

RA-RA*R£FLNO 

R>»RA 

30  CONTINUE 
IFLAG-1 
GO  TO  70 
40  CONTINUE 

REFLNO-REFL (ANGLN ) 

IF  (N.GE.2)  GO  TO  50 
RB-REFLNO 
R*RB 
GO  TO  60 
30  CONTINUE 

RB-RB*REFLN0 

60  CONTINUE 
IFLAG-0 
70  CONTINUE 


C  CALCULATE  PARAMETERS  Z  AND  PRES 
C 

Z-DCMPLX <  0 . ODO . UN*DISTX*DCOS ( ANGLN >  > 
PRES«CDEXP<Z)*<-1>**INTUN+1. )/2. >* 

C  ( 1 . ODO+1 . OPO/REFLNO)*R 
C 

C  CALCULATE  THE  PRESSURE  P  WHICH  IS  A  SUM  OF  PRES 
C 

P-P+PRES 

C 

C~5ET"TD~THE"HEXT-PAIR  OF'IMAGE 
C 

N-N+l 
GO  TO  10 
80  CONTINUE 
C 

C  CALCULATE  PA  AND  PH. AND  RETURN  BACK  TO  THE  MAIN  ROUTINE 
C 

PAMP-CDABS(P) 

PREAL-P 

P-P-PREAL 

PP«P*DCMPLX  <  0 . ODO . 1 . ODO ) 

PIHAG»DSIGN(CDABS  <P>  »-PP> 

PHAS-DAT AN2 ( P IMAG . PREAL ) 

RETURN 

END 


non  nnn  ooono  noon  non  non  noon 


*«*  SUB  PROGRAM  *** 

( (PROGRAM) ) 

FUNCTION  REFL  ( ANGLN) 

DEFINE  CHARACTERS  OF  VARIABLES.  PARAMETERS  AND  CONSTANTS 

IMPLICIT  REAL*8 < A-H » 0-Z ) 

COMPLEX* 1 6  REFL. PS A I 
COMMON  C21.RC21 

CALCULATES  PARAMETER  CHECK 

CHECK-1. 0D0-C21**2SDC0S< ANGLN )**2 

IDENTIFY  THE  GRAZING  ANGLE  WHICH  IS  LESS  THAN  THE  CRITICAL 
ANGLE  OR  NOT 

IF  (CHECK. GT. 0.000)  GO  TO  10 

CALCULATE  PARAMETER  PSAI  FOR  THE  REFLECTION  COEFFCIENT  OF  EACH 
BOUNCE.  THERE  ARE  TWO  WAYS  WHICH  DEPEND  ON  THE  IDENTIFICATION 
OF  A  CHECK 

PSAI-  DCMPLX(O.ODO.-OSQRT ( -CHECK >/DS IN ( ANGLN) ) 

GO  TO  20 
10  CONTINUE 

PSAI-  DCMPLX ( DSQRT ( CHECK > /DSIN( ANGLN > . 0 . ODO > 

20  CONTINUE 

CALCULATE  REFL 

REFL-(RC21-PSAI )/(RC21+PSAI > 

RETURN  BACK  TO  THE  SUBROUTINE  PRES 

RETURN 

END 
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